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Abstract 

We propose a class of auxetic three-dimensional lattice structures. The elastic mi¬ 
crostructure can be designed in order to have o mni -directional Poisson’s ratio arbitrarily 
close to the stability limit —1. The cubic behavior of the periodic system has been fully 
characterized; the minumum and maximum Poisson’s ratio and the associated principal 
directions are given as a function of the microstructural parameters. 

The initial microstructure is then modihed into a body centered-cubic system that can 
achieve a Poisson’s ratio lower than —1 and that can also behave as an isotropic three- 
dimensional auxetic structure. 
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1 Introduction 

The Poisson’s ratio u, the negative ratio between lateral and longitudinal deformations, 
is an indication of the capacity of a material to resist to different types of deformation. 
For almost incompressible materials like rubber, most liquids and granular solids, external 
mechanical loads must spend much more work to change the volume rather than changing 
their shape. On the contrary re-entrant foams [371ES] and some molecular structures m 
[HIES] can easily change the volume homotetically, but they are hardly deformable in shape. 
For a stable material the Poisson’s ratio ranges between the ‘dilational’ lower bound — 1 
and the ‘incompressible’ upper bound 0.5 in the linear isotropic case, while for anisotropic 
materials constitutive stability defines a domain in a properly defined n—dimensional space, 
as shown in BB and in [IH] for the cubic case. For a perfect ‘dilational’ material a dilation 
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is the only easy mode of deformation [ISl SS], which results in an isotropic material with 
u = — 1 . 

Auxetic materials are materials with negative Poisson’s ratio, expanding (contracting) 
laterally in one or more directions, when stretched (compressed) longitudinally. The term 
auxetic, from the Greek au^rjaic (auxesis: increase, grow), was hrstly proposed by Evans 
[23] • Such non standard materials are increasingly used in sport applications (impact 
protector devices, better conformability for comfort and support, and enhanced energy 
absorption for lighter and/or thinner components) [SD EHl IS], in biomedicine (efficient 
design of dilators and stents, arterial prostheses and bandages) [I?] [2Sl SIl El HD], fabric 
and textile PIIIISSI and military applications |39j . The strong interest for engineering 
applications is related to enhanced properties, such as shear resistance [SSI 133], indentation 
resistance (2D11SI, synclastic curvature [23], crashworthiness [SE], pore-size tunability PEI 
and sound absorption and vibration damping [S3l 11211S21 ISO] • Recently, negative Poisson’s 
ratio materials have been also associated with metamaterials |2ZlEnilSZ|; concerning such 
parallelism, it is important to correct a general misconception that identifies metamate¬ 
rials to all manmade microstructure media; properly speaking metamaterials concern the 
dynamic behavior in the low-frequency regime IZB 

While Love in his famous treatise [3S| already mentioned materials with negative Pois¬ 
son’s ratio, they started to become popular after the first artificial examples, such as the 
re-entrant three-dimensional structure of PES], the re-entrant polymeric foam of 1371 . and 
the continuum models of [321133] • Extended reviews of existing models can be found in 

[3311221SH SZj. 

A mechanism made of rotating rigid units can exhibit auxetic behaviour as proposed 
by HDj; it is composed of rigid cuboids connected at their edges, which deform through 
relative rotation with respect to each other. The system exhibits negative values for all the 
six on-axis Poisson’s ratios. The use of three-dimensional models to predict or to explain 
auxetic behaviour is used in PPil where a three-dimensional rotating and/or dilating 
tetrahedral model is applied to real crystalline auxetic materials such as «-cristobalite and 
«-quartz structures of both silica and germania. Auxetic responses have been demon¬ 
strated possible in many crystals [131 [I3|, and some artificial three-dimensional auxetic 
materials have been proposed especially with the advent of the 3D-printer. A new ma¬ 
terial, named ‘Bucklicrystals’, achieves a three-dimensional auxetic behavior through the 
elastic buckling [Hj. In [IH] it is described a dilational three-dimensional cubic auxetic 
material with an ultimate Poisson’s ratio of v = —1, based on a two-dimensional chiral 
model recently published |3S|; the three-dimensional model is studied numerically and a 
sample is fabricated with a 3D-printer. The development of the modern 3D-printers has 
allowed the creation of metamaterials starting from theoretical models with unit cells in 
the micro-meter range: in ra the authors fabricated and characterized a truly three- 
dimensional crystalline mechanical metamaterial with unit cells in the micro-meter range 
and with adjustable Poisson’s ratios, including negative values. They modified appropri¬ 
ately the DLW technology (Direct Laser Writing) to create three-dimensional nano- and 
micro-structures height only some tens of micro-meters, to get larger structures, and they 
apply this new approach to a model previously introduced in [28] . 

Recently we have proposed a class two-dimensional lattice models having the Poisson’s 
ratio arbitrarily close to —1 [TB|; in this work we rediscovered a cubic model firstly given 
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in isHiisn], where an advanced topology optimization procedure was used in order to de¬ 
sign the optimal structure. In [T^ we also generalized the cubic structure to an isotropic 
structures with omni-directional negative Poisson’s ratio which remains arbitrarily close 
to the stability limit —1. In the model internal hinges avoid concentration of deformations 
that can easily lead to plastic deformations. In addition, the superposition of two systems 
rotating in opposite directions balances any internal couple canceling any chirality. Mi- 
crostructured media with omni-directional negative Poisson’s ratio have also been proposed 
in (65| for plane structures. Here, we extend such plane structures to a three-dimensional 
cubic lattice as in and we analyse in detail the effective behavior considering elastic 
elements and additional internal springs leading to overall stability. The initial system has 
a simple cubic unit cell with strong directional dependance of the effective constitutive 
properties. A modified body-centered cubic microstructure is designed in order to have 
omni-directional negative Poisson’s ratio. 

The paper is organized as follows. In Section we show the simple cubic auxetic lattice 
and we analyse the kinematics of the unimode structure. The effective properties are 
determined in Section 2T and 2^, where we consider hrst the behavior along the principal 
directions of the cubic system and then we search for the extremal values of the Young’s 
modulus and the Poisson’s ratio. The structure shows strong directional dependance and 
a modified body centered cubic structure, stiffened with respect to shear deformations, is 
presented in Section]^ Such modihed structure can be tailored to give isotropic behavior 
or Poisson’s ratio less than —1 along particular directions. Final considerations conclude 
the paper. 


2 The auxetic lattice 

The typical cell of the lattice is presented in Fig. [T^. It is composed of slender cross¬ 
shaped elements, represented in red and blue, joined by a hinge at the central point, as 
in Fig. [T)d. Additional trusses, represented in green, connect the cross-shaped elements 
at the end points. Introducing the terminology of crystal structure, the unit cell is a 
cubic P (cP) Bravais or space lattice jM] with cubic symmetry. In Fig. we indicate 
the geometrical parameters characterizing the structure: p is the length of the arms of 
the cross-shaped elements and 7 their inclination with respect to the truss elements. The 
longitudinal truss elements have cross section At, Young’s modulus Ef and longitudinal 
stiffness kt = EtAt/{2p cos j). The arms of the cross-shaped elements are considered as 
Euler beams with cross section A^ second-moment of inertia Jc and Young’s modulus Ec- 
We tesselate the space starting from the unit cell in Fig. [^. The corresponding periodic 
structure, given in Fig. Il>. has a Bravais periodic lattice m consisting of points 

R = niti-I-n2t2-f nats, (1) 

where ni,n 2 ,n^ are integers and 

/ 2 psin 7 \ ( ^ \ f ^ \ 

ti = 0 , ta = 2 psin 7 , tg = 0 , ( 2 ) 

V 0 / V 0 y \2psm^J 
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Figure 1: (a) Three-dimensional cell of the cubic lattice, (b) Particular of a face of the 
typical cell of the cubic lattice. The numbers identify the node in part (a) and (b). 



Figure 2: (a) Unit cell implemented in Comsol Multiphysics®, (b) Three-dimensional 
cubic lattice. 


are the primitive vectors spanning the lattice. 

In the absence of trusses and for rigid cross-shaped elements we refer to m Section 
3(a) and to [12] for the description of the kinematics of each couple of elements joined 
at the central points. In particular, the only possible mode of deformation is homothetic 
expansion or contraction. Following the systematic analysis for hnite deformation given in 
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( 3 ) 


[121 mi SSI ss] we show that the lattice is a unimode material. Let 

T = [ ti t 2 ta ] = 2psin7l 

be the ‘lattice matrix’, where I is the identity matrix. Dnring the deformation the primitive 
vectors nndergo an affine transformation and the matrix T describes a motion starting at 
t = to, with 7(to) = To- At time t the deformation gradient is given by 

F((,(„) = [T(f)][T(f„)l-V (4) 

Then, the corresponding Cauchy-Green deformation tensor is the path 

C(i,(„) = [T((„)]-’’[T(()F[T(t)||T(i„)l-‘ = I. (5) 

As for the planar square lattice in ra the only possible path C{t,to) lie on a one¬ 
dimensional curve and we can conclude that the lattice structure is unimode. 


2.1 Effective properties 

The lattice structure has cubic symmetry and its macroscopic behavior is described by 
three independent elastic moduli. The elements of the lattice are sufficiently slender, 
so that classical structural theories can be conveniently applied to analyse the effective 
response. The structure has been studied numerically with the finite element code, Comsol 
Multiphysics®. We have considered the following material and geometrical parameters: 
steel with Young’s modulus Ec = 200000 MPa, cross-shaped elements with length of the 
arms p = 10 mm and circular cross section having radius r = 0, 25 mm, so that the area 
is Ac = 0,196 mm^. The truss elements have longitudinal stiffness kL and we introduce 
the non-dimensional stiffness ratio parameter a = k^p/{EcAc), between the longitudinal 
stiffness of the trusses and of the arms of the cross-shaped elements. 

The single unit cell has been implemented, subjected to periodic boundary conditions. 
In particular, displacements u = {ui,U 2 ,u^) along the principal directions ti, 62 , 63 of the 
lattice have been applied on the boundary and the corresponding reaction forces on the lat¬ 
tice have been obtained. Displacements and reaction forces are associated to macroscopic 
strains e and macroscopic stresses a, respectively. The effective Poisson’s ratio 


_ _ ^’' 22^11 ~ ^’' 11^22 _ 

(dll ~ ^ 22)^33 + (dll -|- d22)(dii — £22) ’ 

is given in Fig. as a function of the stiffness ratio a. 0 It is shown that the lattice has 
auxetic behavior in the whole range of a. In the limit when a —7- 0 the Poisson’s ratio u* 
is arbitrarily close to the stability limit — 1 . 

The Young’s modulus 


_ {^11 — ^22)[^ll{^ll — 2622 + ^33) ~ Cr22(~2£ii -|- £22 + ^33)] 
(dll ~ d22)£33 + (dll + d22)(£ii — £22) 


(7) 


^Note that equivalent definitions of v*, E* and ji*, are given by permutation of the indexes 1, 2, 3 in Eqs. 
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Figure 3: Poisson’s ratio ly* as a function of the non-dimensional stiffness ratio a — 

kLp/{EcAc). 


and the bulk and the shear moduli 


K* 


Tr(tT) 

3Tr(£)’ 




<712 


( 8 ) 


are shown in Fig. 

We note that in the limit a —> 0 the Young’s and the bulk moduli tend to zero while 
the shear modulus remains finite (and positive). This is clearly fully consistent with the 
linear theory of elasticity for cubic materials and also observed asymptotically for the 
plane model in [H]. In particular, the elasticity tensor in case of cubic symmetry can be 
expressed as 

C = 3K*I + 2fi*K + 2p*K, (9) 

where p* is a second shear modulus for shear inclined at 45° with respect to two principal 
directions of the lattice [SOI Eg. The fourth-order projection tensors in (|^ are J = (1/3)1® 
I, with I second-order identity tensor, K = T — J and K = I — T. The tensors I and T have 
components 


hjki - ^{SikSji + SuSjk), Tijki - SijSjkSki, (^,T k, I not summed). (10) 

The Poisson’s ratio is then given by 


u* = 


3K* - 2jl* 
2{3K* + fi*y 


( 11 ) 
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Figure 4: (a) Effective Young’s modulus E* and effective bulk modulus K*. (b) Effective 
shear modulus ji*. Elastic moduli are given as a function of the non-dimensional stiffness 
ratio a. 


confirming that v* ^ —1 when K* 0 . 

From the results in Fig. we also note that the shear modulus ji* is independent of the 
ratio a, this is due to the fact that during shear deformation the longitudinal truss elements 
are not deformed and the effective stiffness of the system is associated to longitudinal 
deformations of the arms of the cross-shaped elements, a deformation mechanism already 
observed in the plane model [ TR1| . 


2.2 Directional dependance of effective properties 

In the previous Section the effective properties have been shown in the principal system of 
the lattice. Here we detail the directional dependance of the elastic moduli, in particular 
we focus on the Poisson’s ratio and the Young’s modulus. It is convenient to consider the 
compliance tensor S = which, for cubic material, in addition to the major and minor 
symmetries ^kiij ^nd Sjiki ^ijik} has 

'S'llll = S2222 = S3333, S'1122 == <S '2233 = >S'3311, *S'i212 = <S '2323 = 'S'3131, ( 12 ) 

with the constrained of positive definiteness |61| 

'S'llll — 'S '1122 > 0, 'S'llll + 2 S'ii22 > 0, 'S '1212 > 0. (13) 


In a rotated frame of reference, where the rotation is described by the proper orthogonal 
tensor Q^, the new components of the compliance tensor are 


^ijkl QipQjqQkrQlsS'^ 


-'pqrs • 


(14) 
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The elastic compliance and S[i 22 

‘S'nii = nifijnkniSijki = 5'iiii(ni + + nj) + {2Sn22 + Si 2 i 2 ){nlnl + nlnl + n^rig), (15) 


S[i 22 = nirijmkmiSijki = Snn{nlml + + nlmf) + Sii22{rn\ril + + m\nl + 

mini + +mlnl + mini) + Si2i2{nin2mim2 + nin^mim^ + n2n2,m2m-i), (16) 

where n and m are orthogonal unit vectors. 

The dependance of the Young’s modulus on the direction, E* = £'*(n), can be obtained 
from the relation 

^ (17) 


£'*(n) 


which, making use of Eq. (15) and the identity 

4i 4i 4 1 o/22i 22i 22\ 

+ ^2 + % = 1 + ^2^3)5 

gives 


2_2\1 -1 


£'*(n) — — {2Siiii — 2 S 1122 — Si 2 i 2 ){Ti\n\ + n\n\ + Ti\n\)^ 


(18) 


(19) 


The directional dependance of the Poisson’s ratio z^*(n, m) is obtained from the relation 




* / \ i 

V (n,m) = - — 

O 


1122 


( 20 ) 


1111 


We rewrite S[i22 in Eq. ( 16 ) making use of the following identities: 

1 / 22 i 22 i 2 2 \ / 22 i 22 i 22 i 22 i 22 i 2 2 \ 

1 — [n^mi + 1121712 + n^m^) = (iii;^ii 2 + m2ni + min^ + m^n^ + 1112113 + 11 I 3 II 2 ), 
{n\m\ + nlml + nlml) = — 2 (iiiii 2 iiiiiri 2 + nin^mim^ + 11211311121113 ), 

so that we have the simplified expression 

'5*1122 ~ <51122 + (<5llll — *S'ii 22 ~ <5i212/2)(i 1^111^ + U^Ul^ + II 3 III 3 ). 

Therefore, the directional dependance of the Poisson’s ratio i 7 *(n, m) is 

* . . _ >5(^22 _ <51122 + (<5llll — <S'ii 22 ~ < 5 l 212 / 2 ) (u^lU^ + U^Ul^ + II 3 III 3 ) 

p (n, mj — 

O- 


1111 


S'liii — 2[Siiii — S 1122 — Si2i2/2){n\n^ + + n^n^) 


( 21 ) 


( 22 ) 


(23) 


2.2.1 Effective Young’s modulus 

The Young’s modulus E* is shown in the polar plots in Fig. [^as a function of the direction 
n. For low stiffness ratio a, which is equal to 0.000515 in Fig. the structure has strong 
anisotropy highlighted by the strong variation of F*, namely the ratio E^^^ ^ 1. 

Within the cubic behavior, the anisotropy can be quantified by the dimensionless Zener 
anisotropy factor HIEEI 

( 24 ) 


4S- 


1212 
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Figure 5: (a) Polar representation of the effective Young’s modulus £'*(n). (b) Polar 

diagram of the effective Young’s modulus E* in the [laO] plane. Results are given for 


tt = 0.000515. 



Figure 6: Effective Young’s modulus E*, along [100], [110] and [111] directions, as a function 
of the non-dimensional stiffness ratio a. (a) Effective Young’s modulus for a G (0,0.5). 
(b) Effective Young’s modulus for a G (0,0.001). 


which is equal to 1 for isotropic materials, while (3cub = 0.325 • 10“^ in Fig. 

The extreme values of the Young’s modulus are EJ'qq, and where the subscripts 
indicate the crystallographic directions. They are represented in Fig. as a function of 
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the stiffness ratio a. In particular, the global extrema are 


f^cub ^ 1 • 


l^cub ^ 1 • 


_ _ 

^max ^100 Q ’ 

^1111 

^_ 

^min "^111 n i o C i /I C ’ 

^1111 + ^^1122 + 4^1212 

E*=E*,= ^ 


'^max - ^111 ~ c I oc I /I c ’ 
^1111 + ^^1122 + 4^1212 

^ 

^min -^100 q "> 

^1111 


(25) 


as 


confirmed by the results in Fig. where Pcub < 0.325 • 10 


2.2.2 Effective Poisson’s ratio 


While for linear isotropic stable materials Poisson’s ratio is strictly bounded between — 1 
and 1/2, such bounds do not exist for anisotropic solids, even for those ‘closest’ to isotropy 
such as the cubic material. In |61j it is demonstrated that arbitrarily large positive and 
negative values of Poisson’s ratio could occur in solids with cubic material symmetry along 
specific directions. The key requirement is that is very large relative to other direc¬ 
tions, and, as a consequence, the Poisson’s ratio for stretch close to but not coincident with 
the [111] direction can be large, positive or negative. [SI] also replaced the conventional 
belief that the extreme values of u* are associated with stretch along the face diagonal 
[110] direction (see, for example |14j). 

The Poisson’s ratio i^*(n, m) depends both on the direction of the applied uniaxial load 
n and on the transverse direction m, with m • n = 0. Indeed, to hnd the extreme values, 
it is necessary to consider general n and m; for the purpose of illustration the dependance 
of p* on the direction m for three different n is shown in Fig. [^ where we have introduced 
the polar angle (5 = arccos(n 3 ) and the azimuthal angle (p = arctan(n 2 /ni). 

For the analysis of the maximum and minimum Poisson’s ratio we refer to |61[ 148] . In 
particular, if n is supposed to be given, m depends on just one parameter and so the m 
that gives the extremes of m) is: 


cosd 


^2 

sin0 

-nin^ 

-rii 

+ ! -^ 

-77.2713 

_ 0 

\/i 

_1 - nl_ 


% 7^ 


( 26 ) 


where 6 is arbitrary. 


Inserting (26) into 


leads to the condition of stationarity of 


tan 29 


nin2n3{nl - nfj 
n?ni(l + ni)-ni(l-ni)2’ 


(27) 


where, if 0 is a solution of (27), so is 6 + w/2. 

Following symmetry considerations the analysis can be reduced to the irreducible spher¬ 
ical sector delimited by directions n = 100, 110, 111, which corresponds to l/48th of the 
unit sphere surface as shown in Fig. [^ where A, B and C correspond to n = 100,110 and 
111, respectively. In jlSj it is shown that there are no stationary values of z^*(n, m) inside 
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Figure 7: Polar diagrams of the Poisson’s ratio u* as a function of the transverse direction 
m. Results are given for three different values of n described by the polar angle /3 and 
the azimuthal angle 0, (/5, 0) = (7r/2,7r/6), (tt/S, Svr/S), (vr/d, vr/T) and for stiffness ratio 
a = 0.000515. 



Figure 8: The irreducible spherical sector defined by the vertices corresponding to, n = 
100,110 and 111. Note that the Poisson’s ratios on the edges C—D and C—A are equal. 
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the spherical sector of Fig. but only on its boundary. Hence, we follow the local mini¬ 
mum and maximum Poisson’s ratio when n moves along the path A—B — C—A visualized 
in green in Fig. while the critical directions of m making u* stationary are obtained 
from Eq. (27). 

We show in Fig. [^the local minima and maxima along the path A—B—C—A among which 
it is possible to hnd the global minimum and maximum. Results are given for a = 0.0005 
and 0.05. In both cases the global minimum is at the point A, where n = 100, and it is 
independent on the orthogonal direction m = Oab, as can be easily checked from Eq. (23). 
The global maximum is at the point B, with n = 110 and m = lIO, corresponding to 
0 = 0 in Eq. ([2^. 



Figure 9: Maximum and minimum Poisson’s ratio C for different position of n along the 
path A — B — C — A shown in Fig. Results are given for a = 0.0005 (black lines) 
and a = 0.05 (grey lines). For a = 0.0005 the extreme values are = —0.967 and 
^max — 0.925. For a — 0.05 the extreme values are = —0.232 and i^^ax — 0.998. 

In Fig. I^we also note the high degree of anisotropy of the lattice, the Poisson’s ratio 
increases quite rapidly moving from the direction n = 100. In the direction n = IpO, with 
0 < p < 1, along A—B, the behaviour is always auxetic along specihc transversal directions, 
while in the proximity of n = 111 (point C) the behaviour is not auxetic. 

For a generic stiffness ratio a, the five possible candidates for the global minimum and 
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maximum of z/*(n, m) are the following |48j : 


i/J'io = i/=^(110,110) = - 


2S 


1122 


<5iiii 


5'ii22 + 25'i212 

>S'll22 ~ 25 'i212 


2<S'i212 
>S'ii22 
'S 




Z/2 = l^*(llp2,P2P22) = 


1 { 2»S'iiii + 45'ii22 ~ 4)S'i 212 
'S'llll + 25 'ii22 + 45'i212 
3*S'i212 


+ 


‘S'llll + <S'll22 

12*=l/*(100,-) = - 

1‘S'llll + 25'ii22 + 45'i212| 


1111 


+ 


'S'llll + 2S'ii22 + 45'i212 I'Sllll + 2S'ii22 + 4S'i212| ’ 


where 


Pi = 


'z2* + l/2 
ul - 1/4’ 


P2 = 


1^2*-3/4’ 


Cl — 6[S'im + 3S'iiii5'ii22 + 2 'S^i 22 — 2(5'iiii + 45'ii22)'Si212], 
C 2 = q[3'Siiii + 1 I'Sllll S 1 I 22 + 1051^22 + 2(3,Siiii + 4 S'ii22)5'i212]- 

O 


(28) 


(29) 


The entirety of possible cubic materials with positive definite strain energy are limited by 
the constraints 


— 1 < uliQ < 1, and 


2 ^^ ^110) ^ ^001 


< 1 - ^no- 


(30) 


These two relations define the interior of a triangle in the (i^ooi; ^iio) The interior of 

the triangle represents the entirety of possible cubic materials with positive definite strain 
energy. 


z/* • 

min 

n 

m 

condition 1 

condition 2 

V 

0 

no 

001 

^001 > 0 

n*io > ^001 

1/2 < 

no 

no 

nio ^ —1/2 

^llo ^ ^001 

1 

A 

100 

arbitrary 

^001 < 0 

Klo > ^001 

— (X) < 

iipi 

ilo 

Ao < - 1/2 

n*io < ^001 


Table 1: Global minimum of Poisson’s ratio for cubic material based on the values 
z /^01 and 


Tables and 1^ list the values of the global minimum and the global maximum 
respectively, for all possible combinations of elastic parameters. Each row of the Tables 
and identifies a sector in the triangles shown in Fig. 10 and defined by Eq. (30). These 
sectors define the global extrema for every point within the triangles, so to identify them 
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z/* 

max 

n 

m 

condition 1 

condition 2 

n* < -1/2 

iipi 

llo 

Ao < -1/2 

^1*10 > ^001 

A 

o 

100 

arbitrary 

^001 ^ 0 

1 — 1 

V 

O 

^110 ^ ^ 

no 

no 

nio ^ “3/2 

^llo ^ ^001 

^001 < 3/2 

no 

001 

0 < ^001 < 3/2 

^*110 < Koi 

^2 < oo 

np 2 

P2P2^ 

^001 > 3/2 



Table 2: Global maximum u* of Poisson’s ratio for cubic material 


we need to know the two quantities ^001 and I'liQi from which it is possible to go back to 
the values sought. In Fig. 10 we show the values of ^001 and for our lattice, each dot 
representing a different value of a. For a —>■ 0, (uqq^, ^ (—1, —1), where K*/n* 0, 
while for a —?• oo, (^001) ^llo) (0,1), where K*/ji* oo. For finite values of the stiffness 
ratio a the effective behavior is such that the corresponding points in the (^001) ^llo) plane 
are close to the constitutive stability limit ^001 + i(l — = 0, where jl* oo, but 

inside the triangle, as detailed in the inset of Fig. Hfi. The minimum Poisson’s ratio for 
different a is always given by Uq, while the maximum is for a < 0.326 • 10“^ and for 
a > 0.326 • 10-^ 



Figure 10: Global extrema of the Poisson’s ratio z^*(n, m) for a cubic material in the 
(^001) ^llo) plane. The dots represent the values of the cubic lattice for different values of 
the stiffness ratio a. (a) Global minimum (b) Global maximum 

In Fig. [D we show the global minimum and maximum of the Poisson’s ratio as a 
function of the dimensionless parameter a. The behaviour of the lattice is always auxetic 
when loaded along one of the principal directions n = (100), (010), (001) of the lattice, as 
shown by the minimum Uq and also in Fig. Non uniform omni-directional negative 
Poisson’s ratio is achieved in the narrow interval 0 < a < 10“^, indeed decreasing the 
stiffness ratio a to 0 the global minimum and maximum of n* tend to —1, which is a 
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(a) (b) 

Figure 11 : Global minimum and maximum Poisson’s ratio z/* as a function of the dimen¬ 
sionless stiffness ratio a. (a) a G (0,0.5). (b) a G (0,0.01). 


stability limit. Increasing a the global maximum tends rapidly to 1“ and the global 
minimum tends to 0 “. 

In conclusion, the cubic lattice presents omni-directional negative Poisson’s ratio, which 
varies with n and m, only in a limited range of a; in general, for every a, the behavior 
is auxetic when the unidirectional stress is applied along one of the principal direction of 
the lattice. The stability limit v* ^ —1 can be achieved when a ^ 0, where, in turn, the 
structure becomes isotropic. 


3 Isotropic auxetic lattice 


In the previous section we presented a cubic lattice, where the Poisson’s ratio is negative 
only with respect to some directions n and m, as shown in Fig. |9pTj Here, we modify 


the previous micro-structure and we present a type of three-dimensional lattice that can 
exhibit isotropic negative Poisson’s ratio. In Fig. [^we show the modified lattice, which 
is obtained introducing diagonal beams depicted in magenta and adding a lattice point 
in the center of the unit cell to the eight corner points. Each diagonal beam is mutually 
constrained to have the same displacement at the central point where a hinge is introduced; 
in addition each beam is also constrained by internal hinges at its ends. Doing the same 
similarity with the crystal system used previously, such a lattice structure is a body- 
centered cubic system (cl). We indicate with 77 the ratio between the cross-sectional areas 
of the just introduced diagonal beams and that of the arms of the cross-shaped elements 
W. 


The effective constitutive parameters are analysed as a function of the non-dimensional 
stiffness ratios a and 77 . The directional dependance of the Young’s modulus E* is shown 
in Fig. [ 1 ^ Results are given for a = 0.00515 and each plot corresponds to increasing 
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Figure 12: (a) Three-dimensional unit cell of the isotropic lattice, (b) Previously intro¬ 
duced cubic lattice (depicted in blue and red for the cross-shaped elements and in green 
for the truss elements), (c) System of additional diagonal elements (depicted in magenta). 
Numbers identify lattice points. 


values of the ratio 77 . The comparative analysis indicates a transition from ‘weak shear 
resistant’ structures, with Zener’s anisotropy factor (3cub < 1 to ‘dominating shear resistant’ 
structures, with f^cub > 1 . At the transition point the Young’s modulus is a sphere in the 
{xiX 2 X^) space, indicating isotropic effective behavior, corresponding also to (3cub = 1 - 

To identify the numerical value for which the ratio 77 gives isotropic behavior at each 
a, we compare the shear modulus ii^qq = 77 *, see Eq. a. with the shear modulus = 
-^ 100 /( 2(1 + J^q)) that the lattice would have in case of isotropic behavior. An example for 
a = 0.00515 is given in Fig. B*; the intersection point at 77 = 2.386 identifies the value 
for which the lattice has the desired isotropic behavior. Clearly, this corresponds to the 
value of 77 for which E^qq = E^j^q = The monotonic dependance of the cross-section 

ratio 77 on the stiffness ratio a in order to obtain isotropic behaviour is reported in Fig. 
B). At these values the effective behavior is isotropic; therefore the Poisson’s ratio does 
not depend on the directions n and m and, when plotted along the path A — B — C — A 
as in Fig. the local and global maxima and minima coincide and do not depend on the 
position within the path. 

The isotropic effective Poisson’s ratio as a function of the stiffness ratio a is shown in 
Fig. [ 1 ^ In addition to the property of isotropy, we note that the interval in which the 
behavior is omni-directionally auxetic is strongly increased with respect to the cubic case 
in Section where 77 = 0; more precisely the Poisson’s ratio is negative for a < 0.0325. 
In contrast, it is not possible to reach the lower and upper limits = — 1 and v* = 0.5, 
respectively. 

We conclude the analysis of the lattice structure by investigating the behavior in the 
whole range of 77 for different values of the ratio a. The paths in the (^ 001 ) ^llo) plane are 
shown in Fig. 16 a where the global minimum i^*(n, m) in different regions is also indicated. 
In Figs. [T 6| 3 and 16: we show and ^ function of 77 , for different values of a. 

Isotropic behavior is attained on the fine corresponding to the intersection 

points between and Fig- highlighted with a dashed line. 


It is interesting to focus on the behaviour when 77 ^ 1. In such a case the minimum 
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(a) (b) (c) 



(d) (e) (f) 


Figure 13: Directional dependence of the homogenized Young’s modulus E* for six different 
cross-shaped ratio rj and for a = 0.00515. (a) rj = 0.5. (b) ij = 1. (c) 77 = 2.386. (d) 77 = 5. 
(e) rj = 10 . (f) rj = 100 . 


Poisson’s ratio is < —1, tending to the stability limit /r* ^ 00 on the line = — 1 
when 77 ^ 00 . For a = 5.092 the limiting value is = —1.32. The maximum is — 1 < ^^ 0 * < 0 
for i^ooi — 0) ^001 the interval 0 < i^qoi — P 2 > — ^001 — 2- Omni¬ 
directional negative Poisson’s ratio is obtained for a < 10“^, namely ~ while 


Finally, it appears that increasing a the accumulation point for 77 —)■ 00 is (t'ooi) ^iio) ~ 
(1,-1), corresponding to p^^^ = 7 ^*( 110 . 73 ,110) = — (\/3 + l)/2 = —1.366 and p^^x — 
^001 ~ t- 


4 Conclusions 

We have proposed a new family of three dimensional lattices with negative effective Pois¬ 
son’s ratio. The auxetic behavior is given by the topology of the microstructure and, in 
particular, by the contrast between the stiffnesses (longitudinal and flexural) of the el¬ 
ements of the microstructure. For the cubic structure the constitutive stability domain 
is a triangle in the (z^ooi^Klo) ^^d the minimum and maximum Poisson’s ratios are 

achieved for infinite contrasts, when the boundary of the domain is reached asymptotically. 
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Figure 14: (a) Effective shear moduli /r* and = E'^qo/( 2(1 + i/q)) of the lattice as a 
function of the non-dimensional cross-section ratio r]. Results correspond to a = 0.00515. 
(b) Non-dimensional cross-sectional ratio 77 as a function of non-dimensional stiffness ratio 
a in order to have isotropic behaviour. 



Figure 15: Poisson’s ratio ly* as a function of a in case of isotropic behaviour. 
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Figure 16: (a) Global minimum Poisson’s ratio z^*(n,m) for a cubic material in the 

(^001) ^llo) plane, for different values of the cross-section ratio rj. (b)-(c) Global mini¬ 
mum and maximum Poisson’s ratio z/* for 77 G (1,10000) and 77 G (1,10), respectively. The 
dashed line indicates isotropic behavior. 


The minimum is attained for a —>■ 0 and rj/a ^ 00 and maximum for a —>■ cxo. It is shown 
that the lattice can be designed to have Poisson’s ratio —1, a case of perfect dilational 
behavior. The structure can also be modified in order to have Poisson’s ratio less than 
— 1 along some oblique direction with a limiting value of v* = —{\/3 + l)/2. In addition, 
an isotropic structure is proposed which can have omni-directional and uniform negative 
Poisson’s ratio. 

The microstructured medium is scale-independent, so that it is possible to think to 
technological applications ranging from large scale structural systems as in Givil and Aero¬ 
nautical Engineering, where hinges are actually built, to micro- and nanostructures leading 
to metamaterials. 
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